Families of Matter- Waves for Two-Component Bose-Einstein Condensates 
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We produce several families of solutions for two-component nonlinear Schrodinger/Gross- 
Pitaevskii equations. These include domain walls and the first example of an antidark or gray 
soliton in the one component, bound to a bright or dark soliton in the other. Most of these solutions 
are linearly stable in their entire domain of existence. Some of them are relevant to nonlinear optics, 
and all to Bose-Einstein condensates (BECs). In the latter context, we demonstrate robustness of 

£f} , the structures in the presence of parabolic and periodic potentials (corresponding, respectively, to 

f*"^ ■ the magnetic trap and optical lattices in BECs). 

' The recent progress in experimental and theoretical studies of Bose-Einstein condensates (BECs) [1] has made 
matter-wave solitons physically relevant objects. One dimensional (ID) dark [2] and bright [3] solitons have been 
' observed in experiments, and possibilities for the observation of their multidimensional counterparts were predicted 
^| , [4]. Further study of matter-wave solitons is a subject of profound interest, not only from a theoretical perspective, 
but also for applications, as there are possibilities to coherently manipulate such robust structures in matter-wave 
devices, e.g., atom chips, which are analogs of the existing optical ones [5]. On the other hand, many results obtained 
for optical solitons as fundamental nonlinear excitations in optical fibers and waveguides (see, e.g., recent reviews 
. [6,7]) suggest the possibility to search for similar effects in BECs. 
^sO ' A class of physically important generalizations of the nonlinear Schrodinger (NLS) equation for optical media, or 
its BEC counterpart, the Gross-Pitaevskii equation (GP), is based on their multi-component versions. In particular, 
the theoretical work has already gone into studying ground-state solutions [8,9] and small-amplitude excitations [10] 
of the order parameters in multi-component BECs. Additionally, the structure of binary BECs [11], including the 
£f~) " formation of domain walls in the case of immiscible species, has also been studied [11,12]; ID bound dark-dark [13] 
and dark-bright [14] soliton complexes, as well as spatially periodic states [15], were predicted too. Experimental 
results have been reported for mixtures of different spin states of 87 Rb [16] and mixed condensates [17]. Efforts were 
made to create two-component BECs with different atomic species, such as 41 K- 87 Rb [18] and 7 Li- 133 Cs [19]. 

In this work, we report novel solitons in the context of coupled two-component GP equations. These solutions 
correspond to new families of solitons even for the NLS equations per se, hence they are also interesting as nonlinear 
waves in their own right. We start by demonstrating their existence in the context of two coupled NLS equations. 
Some of them arc relevant as new solitons in nonlinear-optical models as well. We also demonstrate that all the 
solutions proposed herein persist in the presence of the magnetic trap and optical lattice (OL), i.e., parabolic and 
sinusoidal potentials [20], which are important ingredients of experimental BEC setups. 

Assuming that the nonlinear interactions are weak relative to the confinement in the transverse dimensions, the 
transverse size of the condensate is much smaller than its lengths. In this case, the BEC is a "cigar-shaped" one, and 
' the GP equations take an effectively ID form [21]: 

l -dT = ~^--d^ + T. a M \ + Vj(x) Uj , (j = i,2) (i) 

J k=l 

where Uj(x,t) are the mean-field wave functions of the two species, t and x are, respectively, measured in units of 
2/t^ix and the transverse harmonic-oscillator length a\± = ^Jh/ (miwi) (rrij and uij± are the mass and transverse 
confining frequency of each species), while Tluji±/2 is the energy unit. The coefficients ajk in Eq. (1), related to the 
three scattering lengths ajk (note that ot\2 — t*2i) through cijk — / tirmi(ajk/ai±)(m,j + rak) / (rrijmk), account for 
collisions between the atoms belonging to the same (djj) and different (ajk, j ^ k) species; they are counterparts of 
the, respectively, self-phase and cross-phase modulation in nonlinear optics. While in optics only specific ratios of the 
nonlinear coefficients are relevant (such as a^faa = 2 or a^jaa = 2/3 [6]), in the BEC context the interactions are 
tunable [9,15], especially because they can be modified by means of the Feshbach resonance (i.e., by magnetic field 
affecting the sign and magnitude of the scattering length of the interatomic collisions) [22]. The Feshbach resonance 
allows one to switch between attractive and repulsive interaction [23], and even to switch it periodically in time, by 
means of an ac magnetic field, which allows to create a self-confined 2D BEC without the magnetic trap [24]. 

In this work, we consider the case with m\ — m,2 = m and an = 022, which corresponds to the most experimentally 



X 



1 



feasible mixture of two different hypcrfine states of the same atom species, or, approximately, to different isotopes of 
the same alkali metal, trapped in the potential including the magnetic trapping and OL components: 

Vi(x) = V 2 (x) = V(x) = (n 2 /2) x 2 + V sin 2 (kx + <P). (2) 

In Eq. (2), fl 2 = 2lo 2 /uj 2 _ (u>i x = u>2 X = w x are the confining frequencies in the axial direction) and Vq (measured in 
units of the recoil energy [20]) set the respective potential strengths, k is the wavenumbcr of the interference pattern 
of the laser beams forming the OL, and is an adjustable phase parameter (</> e {0, 7r/2}). To estimate physical 
parameters, we resort to a mixture of two different spin states of 87 Rb, confined in a trap with the transverse frequency 
u>i± = 183 rad/s, which implies that the length and time units are 2^m and 5.46 ms, respectively. 

We consider a rather general case, in which the inter-atomic interactions in the first species are repulsive [therefore, 
we will use the normalization an = +1 in Eqs. (1)], while in the other species they may be either attractive or 
repulsive. As concerns the interactions between the different species, they are, typically, repulsive. Nevertheless, in 
the case of two different spin states of the same atom species, the Feshbach resonance between such states is possible 
too (experimental studies of the Feshbach resonance in this case are currently in progress [25]), therefore attractive 
inter-species interactions may be relevant, and this case is also considered below. The solutions reported herein, and 
their existence and stability regimes are summarized in Table I. 
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Table 1: Existence and stability of structures in the binary BEC. In the "existence" column, +/— indicates the 
repulsive/ attractive character of the respective inter-atomic interaction which is necessary for the solution to exist. 
The "stability" column indicates the sign of the coefficient a 2 2 (we normalize an = +1, and set a2 2 = ±1) and an 
interval of the values of ai 2 for which the solution is stable. 



In most cases the existence and stability of the solution families is investigated numerically. The numerical method 
was implemented as follows: we first seek stationary solutions by means of Newton iterations which are applied to 
the steady-state equations [iuj — ~iij tXX + Y^t=i a jk\ u k\ 2 Uj + V(x)uj (ji is the chemical potential). Subsequently, 

we perform the linear-stability analysis of the obtained soliton solutions uf\x), setting the perturbed solution to be 



u^\x) + e (bje lut + Cje lu *) , where lo = u> r + iuji is a (generally, complex) perturbation eigenfrcquency. 

Then, the ensuing linear stability problem [26] is solved for the eigenfrequencies and eigenfunctions {bj, c*}. Whenever 
the solution is unstable, we also examine its evolution in direct simulations of the full equations (1), using a fourth- 
order Runge-Kutta time integrator with the time step dt = 0.005 (27.3fis in physical units). To initiate the instability 
development, a uniformly distributed random perturbation of amplitude ~ 10~ 4 was typically added to an unstable 
solution. 

We now examine in detail the solutions shown in Table 1. First, in the absence of the external potential, a family 
of domain walls can be found in an exact form for the special case, ai 2 = 3an = 3a 2 2: 

Uj {x, t) = Ae'^ [1 + {-iy tanh(r?a;)] , (3) 

where the chemical potential is /x = Aa^A 2 , and rj 2 = 2a^A 2 (they follow the pattern of domain-wall solutions found 
long ago in the context of coupled Ginzburg-Landau equations [27]). These solutions exist only if ai 2 > and \i > 0. 
Similar patterns were found in Ref. [12] and other related structures were also predicted to occur in higher dimensions 
[28]. 

We have confirmed the existence and stability of the domain walls by direct numerical simulations (not shown here), 
using numerical continuation to extend them to the case ai 2 ^ 3an, where the analytical solution is not available. 
We have thus found that the domain walls exist and are stable for values of ai 2 down to ai 2 = 1. The case ai 2 = 2 
is relevant to nonlinear optics; stability of the domain wall family for this case was suggested by recent numerical 
results obtained for a similar discrete couplcd-NLS model [29]. Here, we find that these solutions are robust as well 
for other values of ai 2 , and, as will be shown below, also in the presence of the external potential in Eq. (2). 
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FIG. 1. Left: The top panel shows the number of atoms (in normalized units) in each component (circles and stars) of the 
dark-antidark soliton. The middle and bottom left panels show the spatial profiles of this solution for 012 = 0.2 and a\i = 0.5, 
while the right panels show the respective spectral planes (u> r ,u)i) of the corresponding eigenfrequencies, as found from the 
linearized equations (toi — implies linear stability). Right: The top panels show the number of atoms in each component of 
the dark-gray soliton and the instability growth versus a\2- The middle and bottom left panels show the spatial profiles of this 
solution for ai2 = —0.5 and an = —0.99, while the right panels show the respective spectral planes. 

Proceeding to the other solutions in Table 1, we first note that, upon fixing an = +1, we set 022 = ±1 for the 
repulsive and attractive interactions in the second species. Although setting a.22 to be 1 formally limits the generality 
of the results, it has been checked that taking |a22 1 7^ 1 produces similar results to those displayed here. 

Given that we have assumed the first component as being always self-repulsive, we look for solutions starting from 
the uncoupled limit (012 = 0) by taking an initially uniform distribution in this component, u\ = 1. If the second 
component is also repulsive, we take a dark soliton as its initial configuration. If it is self-attractive, a bright soliton 
is initially set in it. 

The case of the self- repulsion and dark soliton in the second component, with the inter-species coefficient being 
repulsive too, a 12 > gives rise to a stationary antidark soliton (i.e., a hump on top of a nonzero flat background) 
in the first component, see the left panel in Fig. 1. It is easy to understand this structure, as the atoms in the first 
component, being repelled by the matter in the second one, concentrate in an effective potential well generated by 
the dark soliton (void) in it. Antidark solitons are well-known to occur when higher-order effects (such as third-order 
dispersion) or a saturable nonlinearity are present in the single-component NLS equations, which is possible in optics 
[7]. In that case, the antidark soliton is usually described by a KdV-type asymptotic equation (for the elevation on 
top of the flat background) , and is not stationary, running at the respective velocity of sound [7] . The two-species 
soliton with the antidark component, presented here, is the first example of a stationary antidark soliton, that we 
are aware of, in a model without higher-order nonlinearities and dispersions. It is also the first prediction of antidark 
solitons in BECs, which suggests that an experimental verification would be of particular interest. 

If, on the contrary, the interaction between the two self-repelling species is attractive (a.12 < 0; right panel in Fig. 1), 
then the void (dark soliton) in the second component effectively repels the matter in the first one, and thus generates 
a dip, i.e. a gray soliton in it (for a detailed description, sec Rcf. [7]). Such solitons exist in the regular NLS equation, 
but there (as well as in other instances of their presence that we are aware of) they travel at a nonzero speed (the 
faster the shallower the dip is), while here the gray solitons are stationary. 

The bound dark-antidark two-component states persist for < 012 < 0.7, while the dark-gray bound-state branch 
continues down to 012 = —1, getting appreciably unstable in the region —1 < a.12 < —0.83. In the latter case, the 
time evolution (not shown here) leads to formation of moving dark-gray soliton complexes. 

Similar results were obtained for the last two branches of the new solutions: if the second component is now self- 
attractive, then the attractive interaction between the BEC species (022 = — 1, a.12 < 0; see left panel in Fig. 2) gives 
rise to bright- antidark solitons (whose existence is explained by the fact that the bright soliton attracts matter in 
the other component). Such a type of two-component solitons was predicted in optics [30], but there it cannot exist 
without the third-order dispersion. In our case, this solution branch terminates at a 12 = — 1, since both states become 
completely flat. 

Lastly, if the bright soliton repels matter in the other component, it naturally induces a dip in it, thus generating a 
bright-gray soliton in the case of 022 = — 1, a.12 > (right panel Fig. 2). It is noteworthy that this bright-gray branch 
is extremely robust; we were able to follow such stable solutions down to ai2 = — 3. 
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FIG. 2. Left: Same as Fig. 1 but for the bright-antidark solitons. The middle and bottom panels are for a\2 = —0.5 and 
CI12 = —0.99, respectively. Right: The branch of coupled bright-gray solitons. The middle and bottom panels are for ai2 = 0.2 
and 0.5, respectively. 




FIG. 3. Left: The bright-antidark branch (the middle and bottom panel are for 012 = —0.5 and —1, respectively) in the 
presence of the potential V(x) = 0.01a: 2 + 0.5 sin 2 (x), see Eq. (2). Right: The same for the domain- wall solutions (the middle 
and bottom panels are for ai2 = 3 and 2, respectively). 

We have also examined these novel solutions in the presence of the OL [i.e., for Vq ^ 0, fi = in Eq. (2)]. We have 
found them to be typically stable when (f> = 0, and unstable when </> = n/2 (as might be expected, since the solution is 
then placed, respectively, at a minimum and maximum of the potential). However, there are exceptions to this rule. 
For instance, dark-antidark bound solitons for ayi = 0.5 are unstable at <fi = 0, in intervals 0.04 < Vo < 0.055 and 
0.09 < Vo < 0.29. This oscillatory instability, involving a quartet of eigenvalues, will be examined in detail elsewhere. 

We have also identified all the solution branches in the presence of the magnetic trap, as well as in the case when 
both the magnetic trap and the OL are present. The branches are extremely robust in the presence of the magnetic 
trap. In particular, Fig. 3 shows the branches for the bright-antidark bound state (left panel) and domain walls 
(right panel), which are always stable in the presence of the combined potential with Q 2 — 0.02 and Vq = 0.5. In 
particular, for a mixture of two different spin states of 87 Rb, the confining frequencies corresponding to this value 
of £1 arc ui x = 18.3 rad/s and u)\± = 183 rad/s in the transverse and axial directions, respectively. Then, the four 
cases shown in Fig. 3 correspond to the mixture containing, in the first and second species, 2 x 10 4 and 6 x 10 3 atoms 
(a 12 = -0.5), 2 x 10 4 and 3 x 10 3 atoms (a 12 = -1), 7.4 x 10 3 in each component (a 12 = 3), and, finally, 7.6 x 10 3 
atoms in each component (a 12 = 2), respectively. 




4 



FIG. 4. Left: The dark-antidark solution branch for <f> = 7r/2. The branch is unstable at arz < 0.4 (an example is displayed 
for ai2 = 0.2, in the middle panel), and then it gets stabilized (see an example for ai2 = 0.86 in the bottom panel). Right: 
Time evolution of the unstable soliton in the case 012 = 0.2 with a random perturbation of an amplitude ~ 10~ 4 added to the 
initial condition. The instability leads to an oscillating dark soliton. The time unit is 5.46 ms, hence the onset of instability 
occurs at « 0.33 s. 

Finally, we notice that the increasing interaction between the components can play a stabilizing role for solutions 
that are unstable in the OL with <f> = ir/2 (Fig. 4). The dark-antidark coupled state serves as an example, being 
unstable for < an < 0.4 and stable for a\2 > 0.4 (left panel). The instability evolution is shown in the right panel 
for ai2 = 0.2: the dark soliton becomes mobile and starts to oscillate in the combined potential. For the physical 
example mentioned above, the numbers of atoms in the two components corresponding to the latter case are 10 4 and 
1.6 x 10 4 , respectively. 

In conclusion, we have presented a number of novel families of composite solutions in the generic two-species BEC 
model. In particular, the first possibility to create an antidark soliton in BECs is predicted. In most cases, the 
compound solitons and domain walls arc very robust, keeping the stability in the presence of the parabolic and 
periodic potentials. In some cases, the new solutions are relevant to optical models as well. It would be of interest to 
look for such states experimentally. 
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